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Abstract. We study the entanglement of a pure state of a composite quantum system consisting of several 
subsystems with d levels each. It can be described by the Renyi-Ingarden-Urbanik entropy S q of a decom¬ 
position of the state in a product basis, minimized over all local unitary transformations. In the case 9 = 0 
this quantity becomes a function of the rank of the tensor representing the state, while in the limit q —> oo 
the entropy becomes related to the overlap with the closest separable state and the geometric measure of 
entanglement. For any bipartite system the entropy Si coincides with the standard entanglement entropy. 
We analyze the distribution of the minimal entropy for random states of three and four-qubit systems. In 
the former case the distributions of 3-tangle is studied and some of its moments are evaluated, while in the 
latter case we analyze the distribution of the hyperdeterminant. The behavior of the maximum overlap of a 
three-qudit system with the closest separable state is also investigated in the asymptotic limit. 


1 Introduction 


After more than twenty years of intensive research entanglement of the pure states of bipartite quantum 
systems is rather well understood |1,2] for two subsystems of an arbitrary dimension d. In this case any pure 
state can be represented in a product basis by a matrix of coefficients of order d, and its standard singular 
values decomposition allows one to reveal entanglement. 

On the other hand, in the case of systems composed of n > 3 subsystems the problem becomes much more 
complicated, as the state is represented by a tensor of size d and n dimensions. Even though several important 
results were obtained, especially in the case of three [3}|5| and four qubits [6J--10], and several measures of 
entanglement in such systems were proposed 11-14 , it is fair to say that the complete understanding of the 


phenomenon of entanglement in multipartite systems is still missing. 

To characterize entanglement of a quantum state of a bi-partite system it is natural to analyze the degree 
of mixing of the reduced density matrix. For instance, making use of von Neumann entropy, one arrives at 
one of the most often used measures: the entropy of entanglement 15 . It is often convenient to apply for 


this purpose also the generalized entropy of Renyi [2] or some other kinds of entropy. 

The aim of the present work is to propose a possible generalization of this quantity for the case of 
multipartite systems, for which a pure state is represented by a tensor. Furthermore, we would like to 
make a connection with the geometric measure of entanglement, which depends on the distance of the state 


considered to the closest separable state 16-18 


Following the papers of Parker and Rijmen 19 and Bravi 20 we suggest to analyze the entropy of 


decomposition of a quantum state in a product basis, sometimes called Ingarden - Urbanik entropy 21 
minimized over all local unitaries. This quantity can be generalized in the sense of Renyi. 
to establish a direct link with the geometric measure of entanglement 


16 18 


Interestingly, 
it is sufficient to consider the 


Renyi - Ingarden - Urbanik (RIU) entropy and send the Renyi parameter q to infinity. 
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Even though the approach advocated here is applicable for arbitrary composite quantum systems, for 
concreteness we concentrate the majority of this work for the case of three and four qubits. It is demonstrated 
that there is no a single pure state for which the minimal RIU entropy is the largest for all values of the 
Renyi parameter q. Investigating the problem for selected values of q we identify certain pure states, which 
are conjectured to maximize this particular measure of entanglement. 

Furthermore, we make also use of a statistical approach to analyze the distribution of the hyperdetermi¬ 
nant and minimal RIU entropy for random quantum pure states. They are distributed with respect to the 
unique unitarily invariant Haar measure on the space of pure quantum states, the complex projective space, 
CP N ~ 1 , where the total dimension of the complex Hilbert space is N = d n . Analyzing systems composed of 
three subsystems of an arbitrary dimension d we obtain a bounds for the geometric measure of entanglement 
for generic states of such a system. 

This paper is organized as follows. In Section 2 we introduce the RIU entropy for a pure state of a 
multipartite system, while in Section 3 different techniques of tensor decomposition are reviewed. In Section 
4 we present results obtained for three-qubit system, while analogous results for 4 qubits are presented in 
Section 5. A more general case of three subsystems consisting of an arbitrary number of levels is discussed 
in Section 6. Computations of the moments of the distribution of 3-tangle and a derivation of the bound 
for the geometric measure of entnaglement are relegated to the Appendix. 


2 Minimal Renyi—Ingarden-Urbanik entropy 


Consider a quantum state describing a system consisting of n subsystems, with d levels each, \ijj) £ Hn = 
T~L® n . Working in an arbitrary product basis one can represent such a state by a n-index tensor, 


W) = E " ' E ^* 1 . 12 ,-,»«l*i) ® 1 * 2 ) ® • • • <8> \i n ) with C iui2 ,... iin £ 

21 — 1 i n — \ 

The standard normalization condition, (ip\ip) = 1, implies that 


(1) 




= 1. 


( 2 ) 


ii =1 


It will be convenient to introduce a multi-index p = (ii, 12 , ■ ■., i n ), where p can be identified with the 
set {!,..., N = n d }, and use a shorter notation^ = \C^\ 2 = | C'i 1 | 2 . Hence p{\ij))) = (pi,...pn) 


,* 2 , 

represents an iV-point probability vector p, which can be characterized by the Renyi entropy S q (p) = 
y2— 1°® ( EE P'h) ' ^ or 9 1 ttu s quantity reduces to the standard Shannon entropy S(p) = — XE=i ^ 1°S?V> 

which in the context of the decomposition of the state \ip), is called the Ingarden-Urbanik entropy [21 


23 


and written <S' IU (|V’)) = S{p (|V’)))- I n will be convenient to use natural logarithms throughout this paper, 
written log 2 « 0.693. 

The product basis |U) <8> \ii) <8> ■ • • <§) | i n ) is determined up to a local unitary transformation, U\ oc = 
V\ 0 V 2 0 ■ ■ ■ <8 V n , where an unitary matrix Vj acts on the j- th subsystem. As the Ingarden-Urbanik 
entropy of the decomposition depends explicitly on the choice of the product basis it is natural to analyze 
the optimal value, minimized over the set of local unitaries 


22,23 


We shall study a more general case of the Renyi entropies S q with Renyi parameter q > 0. For any 
multipartite state \ip) £ 'H® n we define the minimal Renyi-Ingarden-Urbanik (RIU) entropy, 


S. 


RIU 


(tjj) := min S q \p(U\ 0 

Uloc 


;\m, 


( 3 ) 


where the minimum is taken over entire set of local unitary transformations. 

Proposition 1. For any N qudit state, | i/j) £ its minimal RIU entropy is bounded from above, 

E RIU W < logRmax, where 

Umax = d n - nd(d - l)/2. (4) 
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This statement follows directly from the work of Carteret, Higuchi and Sudbery 22 , who showed that 


performing a local unitary transformation U\ oc = U\ ® ® U n , one can always find a product basis such 

that the decomposition ([Tj) contains no more than J? max terms. A suitable choice of n unitary matrices of 
size d allows one to bring d(d — l)/2 elements of the tensor C to zero. Therefore out of all its d n entries at 
least nd{d — l)/2 can be always set to zero. This fact can be also formulated as the following statement 
Proposition 2. For any n qudit state, \ip) £ %® n , represented by a tensor C as in Eq. ([l]), its tensor rank 
R is bounded by R m ax = d n — nd(d — l)/2, where R is the minimal number r of terms in its decomposition 
involving arbitrary coefficients /„, 


cv 


= £A 




(5) 


Observe that in the case of the Renyi entropy of order zero, q = 0, one has 5^ IU (^) = logi? where R 
is the tensor rank of the tensor C{ip). The quantity logi? is known as the Schmidt measure characterizing 
entanglement of the multipartite pure state 24 25 , while the rank R was used to determine probabilistic 
conversion of three qubit pure state [26] . Since the Renyi entropy is in general a non-increasing function of 
the parameter q, one gets the general bound 


Sf U (|V>>)<logW))<logi?„ 


( 6 ) 


For instance, in the case of a three-qubit system, (n = 3 and d = 2), the bound Q gives i? max = 8 — 3 = 5 
in agreement to the five-term standard form of a three-qubit pure state by Acin ei aZ [3j. 

Besides the case a) q = 0, corresponding to the tensor rank of | ip), we shall also consider some other 
particular cases of the definition ©• 

b) q = 1. The minimal IU entropy determines the minimal information gained by environment 

after performing a projective von Neumann measurement of the pure state \ip){ij)\ i n an arbitrary product 
basis [20, 

c) q = 2. The minimal decomposition entropy c h arac t er i zes the maximal purity P 2 = e~ S2 

of the probability vector p^ associated to the outcomes of a projective measurement in a product basis and 
is accessible in a coincidence experiment with two copies of the multipartite state | ip). This quantity was 


used by Parker and Rijrnen 19 to analyze multipartite entanglement in the context of coding theory. 


d) q = oo. In the limiting case q —>■ oo the minimal RIU entropy gives S^ v (\tp)) = — logA max , where 
the largest component of the vector p^ reads A max = max |(V’IXsep)| 2 - The maximum is taken over the set of 
all separable states, |Xse P ) = Uj oc |0- • -0), so A max is a decreasing function of the Fubini-Study distance to 

Dfs = arccos(V%i 


the closest separable state 
entanglement 


18 


16 


17 


, E g (4 0 = 1 - A n 


it), and its function called the geometric measure of 
= l-exp[—S“ u (| m- 


In the case of bipartite systems (n = 2) one can use the standard singular value decomposition of the 
matrix C to show that the probability vector coincides with the vector of Schmidt coefficients A j deter¬ 
mining the Schmidt decomposition !• W = ZU V^IJ>A®I 3)b, where y"f_ 1 \j = 1. This observation 
leads to 

Proposition 3. For any bipartite state, \i/j) £ 7 id <8> 7?<j, its minimal RIU entropy coincides with the 
Renyi entropy of entanglement, so in the case 9=1 one arrives at the standard entropy of entanglement 


S? IU W = E(m = - £ A, log Xj = 5(A). 


(7) 


3=1 


3 Tensor decompositions 

Singular value decomposition (SVD) of the matrix containing the expansion coefficients of a biparite pure 
state plays a crucial role in evaluating the minimal RIU entropy. However, in the case of multipartite systems, 
and quantum states described by higher order tensors, this decomposition is not directly applicable. 
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In this section we consider two generalizations of SVD, the higher order singular value decomposition 
(HO-SVD) 127 and the Parallel Factor model (PARAFAC) [28|. Both decompositions were developed in the 
framework of Principal Component Analysis (29] and have found several applications in signal processing, 
numerical linear algebra, graph analysis, and numerical analysis 30 . Throughout this section we will 


consider such schemes to evaluate some bounds for the minimal RIU entropy and the geometric measure of 
entanglement. 

Let us define a tensor C as a multidimensional or n -way array of numbers, so it can be identified with 
an element of C dld2 '" d ". Accordingly, a matrix and a vector are a 2-way and 1-way tensors, respectively. 
Note that such a space is linear. Given two tensors A,B£ C dld2 '" d " their inner product is inherited form 
the linear space <C dld2 " dn and defined as follows 


(A,B)= ]T Ai 


B, 


l \,%2 5 


( 8 ) 


il^2v)*n 


where the over-bar denotes the complex conjugation. The corresponding induced norm is the Frobenius 


norm 


27 


30 


written ||A|| := (A, A). It implies the Frobenius distance, d(A,B) := ||A — B ||. Note that 
that the coefficients of the state (llj) can be arranged in a tensor C £ C d " and conversly, a given tensor in 
such a space defines a certain pure state \tj}) in the Hilbert space "H?". 


3.1 Higher Order Singular Value Decomposition 

Let C £ (£d 1 d 2 ...d n kg an n _ wa y tensor. We define the fc- th unfolding C^) as the matrix of size dk x 
(■ dk+idk +2 ■ ■ • d n d\d 2 ■ ■ ■ dk- i) that contains the tensor element in the entry with 


j = l + - l )Je, 


with J e = ]Q d m - 


e=i 

i^k 


m—1 

m^k 


The fc-mode product of a matrix £ C dkXdk with C is dehned (element-wise) as 




(fc) 


(9) 


The higher order singular value decomposition allows one to construct a tensor A , called co-tensor 30 , of 
the same dimension than C such that 


A = U (1) <g> U (2) ■ ■ ■ <g> U (n) C, 


( 10 ) 


where each U^ acts according to ^ and any two sub-tensors A ik=p and A ie=q , with p and q fixed, are 
orthonormal 


Ai k=q ) — Sp^q 


r (G 


1 2 


( 11 ) 


(k) 


(k) u 

p are called the fc-mode singular values of C, they are no n-n egative and fulfill a p > a q 


(fc) 


The numbers cr 1 

for all p < q. Such a decomposition is accomplished by taking each U^ in (10) as the matrix of left singular 
vectors of C^ n y Note that finding the SVD of the fc-th unfolding of the coefficients tensor of (|lj) is equivalent 
to diagonalize the reduced density matrix of the fc-th party of the system. Moreover, according to equation 
(101 the co-tensor A of C defines a state |i/’)ho-SVD that is LU equivalent to (JT|) . In this manner, Liu 
et al have recently proposed an entanglement classification based on the study of HO-SVD and the local 
symmetries of the multipartite states [31 . 


3.2 Parallel Factor decomposition 

The idea of expressing a tensor as a sum of rank-one tensors was applied in several contexts. For instance, in 
psychometrics it is known as canonical decomposition (CANDECOMP), while in the brain imaging analysis 
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it is referred to as parallel factor decomposition (PARAFAC) - see 28 and references therein. For other 


applications of this decomposition see also 30 


The PARAFAC is a decomposition of n-way tensor C £ C dld2 ' d " a sum 0 f ran k_ 0 ne tensors, 




u. 


(n) 


i2=k' 


( 12 ) 


fc=l 


where U- !_ k denotes the fc-th column of a matrix of size di x R, with £ = 1 ,,n and A is a vector of size 


R. The symbol o represents the dyadic product of vectors 32 . In terms of the components, the PARAFAC 
decomposition reads 


( 1 ) ( 2 ) 
uU)', ■■■u. 


(n) 


,k' 


Oil ^ ) A k 

k= 1 

In practice, we are interested in an approximation given by sum of r rank-one tensors 


(13) 


C * C PARAFAC = £ A fc UgU o U\% k o • • • 


(14) 


k=l 


Usually a least squares approximation (12) is concern, in which case one has to minimize the quantity 


dp = ||C — ( 7 parafac ||. Since the problem is not linear, usually we obtain an approximate solution only, 
accomplished through the alternating least squares (ALS) algorithm. In this work we use the algorithm 
provided by Nion and De Lathauwer 33 to compute the PARAFAC decomposition of a 3-way complex 
tensor. In the case of r = 1 the state 


/p = 


d 

£■ 

*i=i 


£ Cf£?£\il) ® 1 * 2 ) ® ■ ■ • | in), 


(15) 


is fully separable. Hence, the quantities dp and X p = |(V'|V’p )| 2 are bounds for the geometric measure of 
entanglement Eg and the separable state maximum overlap A max , respectively. 


4 Three qubits 

In this section we analyze pure states of a system consisting of three qubits: n = 3 and d = 2. The exact 
values of the minimal RIU entropy are found for the states W and GHZ for particular values of the Renyi 
parameter q. Generic properties of the minimal RIU entropy are also discussed. Furthermore, some moments 
of the 3-tangle r are computed and the corresponding distribution is analyzed. 

4.1 Minimal decomposition entropy 

In the case of a three-qubit system any pure state | ip) £ Hf 3 = Hg can be represented in the five-terms 
decomposition of Acin et al. |3'|, 

|^) = ai| 000) + o 2 |001) + a 3 |010) + a 4 |100) + o 5 |lll), (16) 

where Yli=i l a *l 2 = anc l f° ur coefficients can be chosen to be real. Observe that selecting ai = as = 1/V2 
and neglecting others, one obtains the state | GHZ) = (|000) +1111))/a/ 2, while setting a 2 = 03 = a 4 = 1/^3 
one has \W) = (|001) + |010) + |100))/V3. As the number of terms in these states cannot be reduced 
by any local unitary transformation their tensor ranks are equal to two and three, respectively, so that 
S^(GHZ) = log 2 and <S ,RIU (1U) = log 3. It is possible to write a state |A 5 ) with all coefficients equal, 
at = 1/y/ft for which S ,RIU (A 5 ) = log5. In the same form, we define |A 4 ) = (|000) + |010) + |001) + 1111))/2 
such that 5 riu (A 4 ) = log 4. 
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For any general value of the Renyi parameter q , we are not aware of any constructive procedure which 
gives the exact value of minimal RIU entropy ([3|. However, for permutation invariant states we can follow 
the general scheme of calculating the maximum overlap with the closest separable state. For such states it 
was first conjectured and later proven that in order to obtain the maximal overlap it is enough to take the 
product state to be a tensor product of the same single-party real state 34-36 . In this spirit, to minimize 


the minimal RIU entropy for permutation invariant states the product of local unitary matrices Ui a 
will be taken as 


U\ oc = U(pr i , with U(p) = 



pe[o,i]. 




(17) 


This task can be done easily in some special cases, i.e. Sj* IU (GHZ) = log 2 for any q and S^ IU (W) = log 3. 
To evaluate the expression ([3| for an arbitrary state, we perform a random walk over the space of unitary 
matrices. Fig.jTJshows the minimal RIU entropy obtained numerically for some particular states and different 
values of q. In the case of the state \W) and |A 5 ) one can compare the value obtained by this procedure 
with the analytical results. 


cRIU 




|GHZ) 

|W> 

|W) Analytic 

|a 5 ) 

IA5) Analytic 

|A 4 ) 

| $ m ax ) 

i<i>r x ) 


Figure 1: (Color online) Minimal RIU entropy obtained numerically for several three qubit states as function 
of the Renyi parameter q. Comparison with analytical results available for states invariant with respect to 
permutation confirms accuracy of the numerical procedure. 


In the following we are concerned with typical properties of the minimum value of RIU entropy. First, 
we compute the value of the Renyi entropy of a random state, then the same quantity is evaluated in 
the corresponding co-tensor and finally we compute the S^ IU by performing the random walk procedure 
described above. As we pointed before, we are mainly interested in the cases q = 1, q = 2 and q = 100. 
The latter value serves as an approximation for the limining case S^, IU . The probability distributions of 
estimations of the minimal RIU entropy for random states are presented in Fig. [2] The mean, second 
moment and standard deviation for these distributions are reported in Table [T] In the case of q = 1 some 


analytical results are available 37 38 


N 


(S 1 ) = *(N + 1)-*(2) = J2 


k—2 


jfc’ 


(18) 


(AS,) 2 = [2tt'(2) ~(N+ 1 )*'(N + 1)], 


where N is the size of the system and 4/ (x) denotes the digannna function and ^'(rr) its derivative. In 
our case we set N = 8 to obtain (Si) = 1.718 and ASi = 0.160, so our numeric calculations are in good 
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agreement with these analytical predictions. We compute the maximum overlap of the state \ip) with the 
closest separable state by performing a random walk in the space of unitary matrices, this quantity is denoted 
as A^x- The PARAFAC decomposition of \ip} yields a bound Ajj^ AFAC f° r this overlap. A comparison 
between these three distributions is presented in Fig. [3} 



<S)j(|'0typical)) 

(| ^typical) HOSVD ) 

S q W (Atypical)) 

q 

<S 9 > 

(s 2 q ) 

A S q 

(s q ) 

(s 2 ) 

A S q 

(s q ) 

(s 2 ) 

<1 

l 

1.717 

2.973 

0.161 

1.132 

1.407 

0.355 

0.907 

0.858 

0.184 

2 

1.534 

2.397 

0.224 

0.812 

0.794 

0.368 

0.650 

0.461 

0.197 

100 

1.125 

1.334 

0.260 

0.488 

0.317 

0.280 

0.383 

0.169 

0.144 


Table 1: Mean value, second moment and standard deviation of the Renyi entropy of the probability vector 
corresponding to a three-qubit random pure state (left), its corresponding co-tensor (center) and the minimal 
RIU entropy for q = 1, 2,100. 


For any value of the Renyi parameter q > 0 one can ask for a state | <F™ ax ), for which its minimal entropy 
S^ lv achieves its maximum value. It is known that in the case of three qubits the maximal entangled 
state with respect the geometric measure is the state | W) [39,40 , and hence the largest value of the RIU 
entropy reads 5'^ TU (IU) = — logA max (IU) = —log4/9 m 0.811. We found numerically states for which 
gRro^rmix) _ i 277, S'| iIU ($™ ax ) = 1.108 attain the maximum values. For q = 1,2 the coefficient vectors 
in the basis © of these maximal states read 


l^max) = o.27| 000) + 0.377| 100) + 0.326|010) 

+0.363|001) + 0.740e-°- 79 ™|lll), 
|$“ ax ) = 0.438|000) + 0.29| 100) + 0.371|010) 

+0.316(001) + 0.698e-°- 8267ri |lll). 


(19) 


Furthermore, we compute also S'^ IU ($“ ax ) for other values of the Renyi parameter q - see Fig. 


0 


4.2 Distribution of 3— tangle 


The residual entanglement or 3-tangle r was introduced by Coffman et al. 42 . It quantifies the genuine 


entanglement of a system of three qubits A, B and C in the following sense. Let Ca,b and Ca,c denote the 
of the density matrices of the pairs of qubits A, B and A, C , respectively. Since the qubits 


concurrences 


41 


BC can be regarded as a single subsystem, we can ask for the entanglement between A and BC. Such a 
quantity will be denoted as Ca,bc and it is equal to 2^/det pa 42 , where pa = trscA4.Bc ' s the reduced 


density matrix of the qubit A when the partial trace respect to B and C has been performed. The following 
inequality holds 


Ca,b + C\,c — Ca,bci ( 20 ) 

and using it we define the 3-tangle as t = C\ BC — C\ B — C\ c . This quantity measures the amount of 
entanglement between the qubit A and the subsystem BC that is not related to the entanglement in the 
pairs A, B and A, C. The 3-tangle r is invariant under the permutation of sub-systems and vanishes on all 
states that are separable under any cut. The 3-tangle for the state (1) with n = 3 and d — 2 is given by 


r(|^))=4|Det 3 (C7)|, 


( 21 ) 


where 


Det 3 (C) 


/ . c «l,3l 


-i2,j2 e kl,il^k2,t2 e i3,k3 e j3^3^h,i2,i3^'jl,j2,j3^--'kl,k2,k3C'ii,g 2 ,£ 3 i 


( 22 ) 


is the Cayley hyperdeterminant of the tensor C. Here , l2 stands for the Levi-Civita tensor of rank 2 and 
the sum is performed over all indexes. The 3-tangle is bounded by 0 and 1. For the | GHZ) state the 
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PCSi) 



P(Sioo) 



P(S 2 ) 




Figure 2: (Color online) Distributions of the Renyi entropy of random probability vector describing a generic 
three-qubit state (□), the Renyi entropy of the corresponding co-tensor (v) and the estimation of minimal 
RIU entropy (o) for (a) q = 1, (b) q = 2 and (c) q = 100. 



Figure 3: (Color online) Comparison between the distributions of ‘S'Ixhj 1 (°), £p = — logA^x AFAC (+) and 
S'lu = — log ALU x ( 0 ) for three-qubit random states. 





























residual entanglement attains its maximum value. Indeed, the pairwise concurrences Ca,b and Ca,c vanish 
and r = Ca,bc = 1 , so this state is referred as a genuinely three-partite entangled state. 

The hyperdeterminant Det3((7) is a homogeneous polynomial function of degree 4, first introduced by 
It is invariant under the action of the group SL( 2,C)® 3 


Caley 43 


44 


This notion of invariance plays 


an important role in the construction of entanglement monotones for pure states of multiqubit systems 45 


Moreover, a classification of multipartite entangled states has been accomplished by analyzing singularities 


of the hyperdeterminant 46 


Now we will discuss properties of typical quantum states of a three-qubit system. Kendon et al. 47 


studied the distribution of 3-tangle for an ensemble of random pure states drawn according to the Haar 
measure. By direct integration of equation (21) with respect the unitary invariant measure they obtained the 
average value, (r) = |. In Fig. [I we show the probability densities of r and r 2 evaluated over a sample of 10 5 
random pure states. Moreover, the first six even moments of r can be calculated by symbolic integration using 
the Beta integral (see Appendix[A]). The distribution of 3-tangle can be approximated by the Beta distribution 


Beta(a,/3;r), using the moment method to estimate the parameters a = (r) 

P=(l~ (t)) ( { $)-{t)2 ~ X ) = ff • Hence we § et 




- 1 


= and 


„ . , „ /31 62 \ 

P b {t) = Beta ^ —, — ;tJ . 


(23) 


By virtue of the chain rule we get an approximated distribution function for r 2 


, „ 1 / 31 62 _\ 

Pb{t) = Vf Beta (l7'l7 ; ^)- (24) 

The above distributions are presented in Fig. [4] together with the estimation obtained by numerical simula¬ 
tions. In order to show the accuracy of our approximation, we show in Table [2] a comparison between the 
first 6 even moments of the distribution (|23[) and the moments of P(r) obtained by symbolic integration. 


k 

fc-th moment of Pb(t) 

(r fe ), k- th moment of P(r) 

1 

1/3 

1/3 

2 

8/55 

8/55 

4 

533/12573 « 0.04239 

128/3003 « 0.04262 

6 

30914/1819783 « 0.01699 

7168/415701 « 0.01724 

8 

112955/13778357 « 0.008198 

98304/11685817 « 0.008412 

10 

1840340/411553533 « 0.004472 

262144/56497545 « 0.00464 

12 

672000151/252556684751 « 0.002661 

4194304/1502700975 «0.002791 


Table 2: Comparison of the first six even moments of r approximated with the equation (23) (center column) 
and computed using the symbolic integration package (right column). 


5 Four qubits 

In this section we compute both the minimal RIU entropy and the hyperdeterminant for several exemplary 
4-qubit states, we also analyze some typical properties of these quantities. 
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P(T) P( T 2 ) 



Figure 4: (Color online) Distributions of the 3-tangle (+) and its square (V) for tliree-qubit pure random 
states. The solid lines are approximations with Beta distributions (231 at left, and (24) at right. 


5.1 Minimal decomposition entropy 

According to Carteret et al [22], a given state |4>) £ Hf 4 = Hie can be written using twelve terms 
|$) = ^2 with c 0 ,1,1,1 = ci, 0 ,i,i = ci,i, 0 ,i = ci,i,i, 0 = 0, 


(25) 


and the other coefficients Cij : k,e are real and non-negative if at most one of the indexes differs from 1. The 
state \GHZ 4 ) is obtained by choosing Co,0,0,0 = c i,1,1,1 = l/v^- Hence, Sq- iv (GHZi) = 2. Moreover, we 
define a state IA42) for which the non-vanishing coefficients in the above expansion are Cij y k,( = 1 /In 
this case S , (( lIU (Ai2) = 12. 

A proper basis for n qubit permutation invariant states is constituted by the so-called Dicke states 

/ \ -1/2 

\D(n,k)} = r\ J2 |1^10^_0). (26) 

' ' permutations n—k 


Note that | W) = |Z?(n, 1)), \W) = \D(n,n— 1)) and | GHZ n ) = (\D(n, 0)) + \D(n, n)))/V2. In computations 
below we will consider the case of n = 4. The calculation of the minimal RIU entropy for permuta¬ 
tion invariant states can be turned into a one-variable optimization by taking Ui oc = U{p )® 4 . We found 
S^ llJ (GHZ 4 ) = log2, for all q and S’{ lIU (D(4,1)) = log4; for the other values of q minimization will be 
accomplished numerically. On the other hand, the optimal decomposition of the state |D(4, 2)) is obtained 
by taking Ui oc = U( 1/2) ® 4 . Accordingly, it is possible to get a compact formula for its minimal RIU entropy 

Sf lv (D{ 4,2)) = log[2 1 - 3 «3- 9 (3 + 3 2? )], (27) 


hence we get 5'( lIU (D(4,3)) = log(8/\/3) by taking the limit q —> 1. On the other hand, when q —>• 00 
one arrives to 5^ IU (D(4,2)) = — log(3/8), which is consistent with the well known result of the overlap of 
|.D(4, 2)) with the closest separable state 34 . Note that by symmetry S'! lIU (D(4,3)) = S'! ?IU (D(4,1)). In 


Fig- E we compare the results obtained in this way with those acquired by performing a random walk over 
the space of unitary matrices for several values of q. 

Consider also some other exemplary 4-qubit states. The hyperdeterminant state 


|HD) = (J1000) + |0100) + |0010) + |0001) + \/2|llll)) , 


(28) 
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maximizes the 4-qubit hyperdeterminant 48 . The minimal RIU entropy for the former state can be com¬ 


puted in the same way as for the Dicke state (27), as it is shown 34-36 that the state (28) presents the 
optimal decomposition minimizing the RIU entropy for any q. Then we obtain 


5f IU (HD) = ^log3, Sf u (HD) 



6 9 \ 
4 + 4 1 ) 


for q^l. 


(29) 


In the limit q^-oowe get (HS) = — log(2/3) which is consistent with the maximum overlap with the 
closest separable states computed according to the Ref. |[34j. On the other hand, the duster states identified 
by Gour and Wallacli [8] 

|Ci) = | (|0000) + |0011) + |1100) - |1111)), 

|C 2 ) = | (|0000) + |0110) + 11001) - |1111>), (30) 

|C 3 ) = i (|0000) + |0101) + 11010) - |1111)), 


are the 4-qubit states that maximize the Renyi a-entropy of partial trace for a > 2. They also found that 
the state 


|l> = vm K 1 + »)(|ooo°) + inn)) + (i - w)( |ooii) + |iioo)) 

+w 2 (|0101) + |1001) + |1010))], with w = exp(2|t), 


(31) 


maximizes the average Tsallis a-entropy of the partial trace for all a > 2. The minimal RIU entropy of the 
former states is shown in Fig. [5] for several values of q. Based on our numerical calculations for the cluster 
states we conjecture that S'! ul +C/ c ) = log 4 with k = 1, 2,3 for any value of the Renyi parameter. 


We also consider the 4-qubit state found by Higuchi and Sudbery 49 


|HS) = i [|0011) + |1100) + w(|0101) + |1010>) + w 2 (|0110) + |1001»], 


(32) 


where w = exp(2f7r/3), which has maximum average von Neumann entropy of partial traces averaged over all 
possible 3 splittings of 4 qubit system into two bipartite systems. Numerical values of 5'^ IU (HS) are shown 
in Fig. [5] We are also interested in the typical properties of the minimum of for four qubits. For an 


qRIU 




|GHZ> 

|D(1,4)> 

|D(2,4)> 

|HD> 

I C*> 

|L) 

|HS) 

|<tj nax > 


Figure 5: (Color online) The minimal RIU entropy S^ lv computed for several four qubit states as function 
of the parameter q. 

ensemble of 10 5 random pure states we compute S q (\i />)), 5' g (|V ; )HOSVD) and S'J lIU (|V’)) taking q = 1, 2,100. 
The corresponding distributions are shown in Fig. [5] and in Table [3] In the first case and with q = lwe can 
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P(Si) 



P(Sioo) 



P(S 2 ) 



-B- 

- s q m) 

-V- 

S q ( I^)hosvd) 

- 0 - 

- sf v m 


Figure 6: (Color online) Distributions of the Renyi entropy of a four-qubit random state probability vector 
(□), the Renyi entropy of the corresponding co-tensor (v) and the minimal RIU entropy (o) for (a) q = 1, 
(b) q = 2 and (c) q = 100. 


compare our results with those provided by the analytic expression (18), (Si) = 2.381 and ASi = 0.124. To 
provide an approximation of S^ IU we analyzed the entropy with the Renyi parameter q = 100. In Fig. [ 7 ] we 
show the distributions of the minimal RIU entropy S^ 1 ^ and S max = — log Alu> where Aj^ ax is the maximum 
overlap with the closest separable state computed through a random walk in the space of unitary matrices. 



'Sjjd'i/’typical)) 

Sq (| ^typical ) HOS VD ) 

5 , r U (IV ; typical)) 

q 

(5,) 

(s 2 ) 

A S q 

(S q ) 

(s 2 ) 

A S q 

(S q ) 

(s 2 ) 

A S q 

1 

2.381 

5.686 

0.124 

2.038 

4.234 

0.283 

1.633 

2.687 

0.145 

2 

2.159 

4.698 

0.190 

1.601 

2.727 

0.403 

1.199 

1.473 

0.192 

100 

1.60 

2.635 

0.254 

1.027 

1.219 

0.405 

0.701 

0.513 

0.145 


Table 3: Mean value, second moment and standard deviation of the Renyi entropy of a four-qubit random 
pure state (left), its corresponding co-tensor (center) and the minimal of the RIU entropy for q = 1, 2,100. 


As in the three qubit case, we are looking for the maximal states of the minimal RIU entropy. For the 
maximal symmetric state with respect to the geometric measure of entanglement |<I> 4 ) = ^/l/3|D(4, 0)) + 
■\/2/3|.D(4, 3)) the RIU entropy yields S^ IU (<b 4 ) = — log(l/3) ~ 1.099 40 . Moreover, our numerical results 
allow us to conjecture that S^^HS) = — log(2/9) « 1.504, is the largest value of the minimal RIU entropy 
with q = 00 among the four qubit states. 

For q = 1 we found numerically a state |4/™ ax ) for which S’^ IU (4'5 nax ) = 1.934 is maximal. In the 
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P(S) 



Figure 7: (Color online) Comparison between the distributions of S'{g I ( ] J (o) and 5 ,LU = — logAj^ x (+) f° r a 
four-qubit system. 


decomposition (25) this state reads 


|$max) = o.630|0000) + 0.28111100) + 0.202|1010) + 0.24|0110) 

+0.232e°' 4947r * 11110) + 0.059| 1001) + 0.282|0101) 

+0.346e-°- 3627ri |1101) + 0.304|0011) + 0.218e°- 626m |1011) 
+0.054e-°- 725,ri |0111) + 0.164e°- 3727 ”|llll). 


(33) 


Further numerical tests support the conjecture that the state |HS), for which S 1 ^ 117 = log 6, maximizes the 
minimal RIU entropy for q > 2. 


5.2 Distribution of the hyperdeterminant |Det 4 1 

In the case of 4-qubit states the hyperdeterminant Det 4 is a polynomial of degree 24. It can be constructed 
following the Schlafli’s procedure [46*. In analogy to the three-tangle we consider the following function of 
the coefficients tensor C of the state | ip) £ ~Hi6, 

T(|V>))= 2 6 3 9 |Det 4 (C)|, (34) 


For a separable four-qubit state one has T(|^ sep )) = 0 but it vanishes also for the states |T>(1,4)), |T>(2,4)), 
|C k ), I HS) and | GHZ). 


Alsina and Latorre found recently [48] that both \HD) and | L) maximize the hyperdeterminant and 
they also discussed relation between these states for which T(\HD)) = T(|T)) = 1. Furthermore, numerical 
simulations indicate that they have the same minimal RIU entropy S)^ 177 . 

We evaluate the quantity (341 over an ensemble of 10 7 random pure states. The mean and the standard 
deviation read (T) = 9.74 x 1CF 4 and (AT) = 2.39 x 10 -3 , respectively, while the corresponding distribution 
is shown in Fig. [8] 


6 Three qudits: Asymptotic case. 


Statistical properties of random tensors and their asymptotic limit became a subject of an intensive research 
50 51 . In this section we study random tenors of dimension three and arbitrary size, which describe 


generic states of a system composed of three qudits. We analyze the bounds for the geometric measure of 
entanglement of such states provided by the tensor decompositions HOSVD and PARAFAC. 
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P(T) 



Figure 8: (Color online) Distribution of absolute value of the hyperdeterminant P(T) (+) for four qubit pure 
random states in a log-log plot. 


Consider a typical state \ip) £ T-L® 3 drawn from the Haar measure. Let A max and A h denote the largest 
component of the probability vectors p and p(IV’)hosvd), respectively, and consider the overlap A p = 
\(ip\ipp)\ 2 , where | ipp) is the state (15). In the cases d = 2,3, we also evaluate the maximum overlap of 
| ip) with the closest separable state by performing a random walk optimization procedure in the space of 
unitary matrices. Fig. [9] (a) shows the average value of four quantities as function of the size of the quit d 
computed over an ensemble of 10 5 random states. The mean value of A max 
the harmonic numbers Hjy, 

N 


can be expressed 52 in terms of 


(A max ) = ^ £ - = H n /N , 
j =i J 


(35) 


where N = d 3 , in our case. Therefore, for a random pure state of three qudits, the average largest overlap 
to a pure state scales as d~ 3 . Performing decompositions for an ensemble of such random tensors we find 
that for large d the largest overlaps optimized by HOSVD and PARAFAC behave as 


(A h ) ~ cT 2 ' 99 , (A P ) ~ d- 1 ' 95 . (36) 

In both cases, the corresponding bound for the geometric measure of entanglement reads (Eq) ~ 1 — (Ak) 
where k = H,P. Fig. [9] (b) shows the geometric measure of entanglement for a three-qudit system as 
function of the qudit size. Moreover, based on the scaling of (Ap) we conjecture that the maximum overlap 
with respect to the closest separable state for a three-qudit system scales as (A max ) ~ d~ 2 for large d. 


7 Conclusions 

We have analyzed the Renyi-Ingarden -Urbanik entropy of pure states of quantum multipartite systems 
minimized over all local unitary operations. For separable states such a quantity is zero irrespectively of the 
value of the Renyi parameter q > 0. In general it is not easy to get analytical results for a given state and 
arbitrary value of the Renyi parameter q. In the special case of permutation invariant states the problem 
becomes easier, as minimization in the space of local unitary operators can be turned into an optimization 
of a one-variable function. 

We computed the minimal RIU entropy of several representative three and four qubit states for various 
values of the Renyi parameter. Some particular states, which maximize the minimal RIU entropy for q = 1 
and q = 2, were identified. Note that this quantity can be considered as a measure of pure states entanglement 
in multipartite systems, and for q —» oo it becomes a function of the geometric measure of entanglement. In 
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Figure 9: (Color online) The mean of the maximum component (left) and the geometric measure of entangle¬ 
ment (right) for random states of a tri-partite system as a function of the qudit size. Bullets correspond to 
numeric simulations: (□) stands for the greatest tensor element A max ; (V) the maximum component of the 
HOSVD co-tensor A#; (o) stands for PARAFAC overlap A p and (+) refers to the overlap with the closest 
separable state maximized by LU. The solid red line (—) is the result (35) with N = d 3 , the solid lines (—) 
and (—) are the best linear fits for HOSVD and PARAFAC, respectively 


the case of three qubits the latter quantity is maximal for the state \W), while for four qubits, it achieves 


maximum for the state (32) of Higuchi and Sudbery 49 


Furthermore, we analyzed the distribution of the minimal RIU entropy for an ensemble of random for some 
selected values of the Renyi parameter, q = 1,2,100. Our numerical simulations demonstrate usefullness 
of the PARAFAC algorithm which allows one the find the principal component of a tensor representing a 
multipartite pure state and estimate its geometric measure of entanglement. 

We studied also the distribution of 3-tangle for random pure states of three qubit system and derived a 
few even moments of this distribution. Numerical results show that this distribution may be approximated 
by a Beta distribution. In the case of four qubits we analyzed the distribution of the absolute value of the 
hyperdeterminant of a random pure state. 

Finally, the behavior of the maximum overlap of a random state of a system consisting of three sub¬ 
systems of size d and the closest separable states was investigated in the asymptotic limit. Although the 
size of the largest component of a random state scales as d~ 3 , numerical results obtained by the PARAFAC 
decomposition of a tensor allow us to conjecture that the minimal overlap scales in this case as d~ 2 . This is 
consistent with an analytic result based on the Marchenko-Pastur asymptotic distribution of singular values 
of random matrices, which gives an upper bound g? _1 . 


8 Acknowledgments 

We are thankful to D. Alsina and J. I. Latorre for fruitful discussions and for sharing their results prior to 
publication. ME is thankful to the Jagiellonian University for kind hospitality and support during his stay 
in Cracow. We acknowledge support of the Mexican National Council for Science and Technology (Conacyt) 
(M.E.) and by the grants number DEC-2011/02/A/ST1/00119 (K.Z) and DEC-2012/04/S/ST6/00400 (ZP) 
financed by the Polish National Science Center. 

A Moments of 3-tangle r 

In the case of even moments of order 2k of 3-tangle r for a three qubit tensor it is easy to note, that it is 
a linear combination of moments of rank 8 k for a random normalized vector. The moments of normalized 
vectors distributed uniformly can be calculated with a use of Beta integral. 
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Let | ijj) be a random unit vector of size d distributed uniformly on a complex sphere, then for a vector p 
of non-negative integers the following expectation value reads 


(Hi\ 2 pi \H 2p2 ■ • ■ \M 2pd ) = m r(Pl J (P2 , +1} • r( f + 1} ■ 

Y(Pi +P2 H- \~Pd + a) 


(37) 


One can also note, that all moments of random vector | ip), which are not in a form presented above, are 

for integrals of monomials over 


53 


equal to zero. This follows, for example, form Collins-Sniady formula 
unitary matrices distributed with Haar measure. 

Calculation of the second moment, (r 2 ) = 8/55, is a simple task, but for higher moments the number of 
term grows so rapidly, that we have used a package for symbolic computations IntU, which yield moments 
presented in Table [2j This package allows for exact calculation of polynomial integrals over the unitary 


group with respect to the Haar measure 54 . 


B Bound for geometric measure of entanglement for tripartite 
states 

The law of Marchenko-Pastur describes asymptotic behavior of singular values of non-hermitian, rectangular 
random matrices. Let X be a random matrix X of size TV x K with entries given by complex random i.i.d. 
normal variables with zero mean and variance 1. We define Y = XX' /txXX' and for c > 1 consider a 
random counting measure on a real line, which counts the number of rescaled eigenvalues of Y which belongs 
to a given set, i.e. 

Mm(A) = ^#{A(cTVY) € A}. (38) 

For a measure defined above, if TV, K —»• oo with additional assumption K/N —» c, there exist a limiting 
distribution /jm —> p given by 


dp{x) = — ^ + ‘ ^ — dx for x £ [a_,a+], (39) 

2tt cx 

with 

a± = (l±v^) 2 . (40) 

The above theorem gives us a behaviour of the largest eigenvalue of a matrix Y. which is Ai(Y) ~ ^ (1+y^) 2 - 
In the case, when TV, K — > oo but K/N —► oo, the theorem does not give us the limiting distribution, but 
form the theorem we will extract, the rate of convergence of the largest eigenvalue. 

Consider, the case when K = TV 2 , then K/N = TV, the direct usage of the Marchenko-Pastur law would 
give us a behaviour of the largest eigenvalue of Y, i.e. Ai(Y) ~ (1 + VTV) 2 ~ . 

Now we will try to use the above asymptotics to bound the geometric measure of entanglement for a 
tripartite random states. 

Consider a random tensor \ip) £ C d ( 8 ) C d <g> C d ='H\® ® FT 3 . The geometric measure of entanglement 

is related to the overlap with the nearest product state, i.e. 

max |(V#)| 2 - (41) 

|0)Gsep 

In the above equation a set sep consist of vectors in a form (f>i (g) c/ 2 ® <f >3 for (f>i £ dii- For distinct i,j , k £ 

{1,2,3}, we denote sep,^ vectors of the form </>i (g) </> 2 for </>i £ T~L, and </> 2 £ Kj < 8 > Hk- We have, that 

sep C sepj| jfc which gives us 

max \{4>\<p)\ 2 < max |(i/|0)| 2 . (42) 

lw 6 se P ld) 6 se Pi|jfc 

The last maximum above is a square of the largest Schmidt coefficient for vector | ip). 
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If we consider, the behaviour of a large random tripartite states the above inequality combined with the 
relations obtained form the Marchenko-Pastur law we obtain 


max |(V#)| 2 < max |(^|</>)| 2 ~ K (43) 

IWSsep |0)esep i!jfc a 
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